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We consider the holographic dark energy (HDE) model in the framework of DGP braneworld 
with Granda-Oliveros infrared (IR) cutoff, L = (aH + /3ff 2 ) -1 ^ 2 . With this choice for IR cutoff, 
we are able to derive evolution of the cosmological parameters such as the equation of state and 
the deceleration parameters, w and q, as the functions of the redshift parameter 2 . As far as we 
know, most previous models of HDE presented in the literatures, do not gives analytically ui = ui(z) 
and q = q(z). We plot the evolution of these parameters versus z and discuss that the results 
are compatible with the recent observations. With suitably choosing the parameters, this model 
can exhibit a transition from deceleration to the acceleration around 2 ss 0.6. Then, we suggest a 
correspondence between the quintessence and tachyon scalar fields and HDE in the framework of 
DGP braneworld. This correspondence allows us to reconstruct the evolution of the scalar fields 
and the scalar potentials. We also investigate stability of the presented model by calculating the 
squared sound speed, i> 2 , whose sign determines the stability of the model. Our study shows that 
i> 2 could be positive provided the parameters of the model are chosen suitably. In particular, for 
a > 1, P > 0, and a < 1, /3 < 0, we have u 2 > 0 during the history of the universe, and so the 
stable dark energy dominated universe can be achieved. This is in contrast to the HDE in standard 
cosmology, which is unstable against background perturbations and so cannot lead to a stable dark 
energy dominated universe. 


I. INTRODUCTION 

A complementary astrophysical data from type la Supernova, Large Scale Structure (LSS) and Cosmic Microwave 
Background (CMB) indicate that our Universe is currently undergoing a phase of accelerating [l|. A component 
which is responsible for this accelerated expansion is usually dubbed ” dark energy” (DE). The simplest candidate for 
DE is the cosmological constant [U which is located in the center both from theoretical and observational evidences. 
However, there are different alternative theories for the dynamical DE scenarios which have been proposed to interpret 
the accelerating universe. One of these models, which has arisen a lot of enthusiasm recently, is HDE. It was shown 
by Cohen et al. [3i] that in quantum field theory a short distance cutoff could be related to a long distance cutoff (IR 
cutoff L) due to the limit sets by black hole formation. If the quantum zero-point energy density is due to a short 
distance cutoff, then the total energy in a region of size L should not exceed the mass of a black hole of the same size, 
namely L 3 pp> < LM 2 V The largest L is the one saturating this inequality, so that we obtain the HDE density as Q 

p D = 3 c 2 M^L~ 2 , (1) 

where M 2 ^ = 8nG is the reduced Planck mass. In general the factor c 2 in holographic energy density can vary with 
time very slowly Q. By slowly varying we mean that (c 2 )/c 2 is upper bounded by the Hubble expansion rate, H , i.e., 

i 

^ < H. (2) 

It was also argued that c 2 depends on the IR length, L [3]. For the case of L = H~ 3 , one can take c 2 approximately 
constant in the late time where DE dominates (fl m < 1/3) JJ. Since in the present work we consider the late time 
cosmology where the DE dominates, we shall assume the factor c 2 to be constant. 

Depending on the IR cutoff, L , many HDE models have been proposed in the literatures. A comprehensive, but not 
complete, list of IR cutoffs which have been used includes the particle horizon radius Rh = cl f* dt/a = a dal (Ha 2 ) 
0,0 , the Hubble horizon L = H - 1 the future event horizon radius Rh = a J t °° dt/a = a / a °° da/(Ha 2 ) |6j, 

the apparent horizon radius L = (H 2 + k/a 2 ) - 1/2 0, the Ricci scalar curvature radius Rcc = (H + H 2 ) - i/2 a 


* mhd@shirazu.ac.ir 
t asheykhi@shirazu.ac.ir 



2 



the so-called Granda-Oliveros (GO) cutoff, which is the formal generalization of Rcc > namely L = ( aH + f3H 2 ) 1 / 2 

0 , the age of our Universe T = / Q a da/(Ha) [15], the conformal age of our Universe V — J 0 dt/a = j ( " da/(a 2 H) 
, the radius of the cosmic null hypersurface [lj|, etc. 

On the other side, in recent years, the theories of large extra dimensions in which the observed universe is realized as a 
brane embedded in a higher dimensional spacetime, have received a lot of interest. In these theories the cosmological 
evolution on the brane is described by an effective Friedmann equation that incorporates non-trivially with the 
effects of the bulk onto the brane. One of the well-know picture in the braneworld scenarios was proposed by Dvali- 
Gabadadze-Porrati (DGP) jl9]. In this model our four-dimensional Universe is a Friedmann-Robertson-Walker (FRW) 
brane embedded in a five-dimensional Minkowskian bulk with infinite size. In this model the recovery of the usual 
gravitational laws on the brane is obtained by adding an Einstein-Hilbert term to the action of the brane computed 
with the brane intrinsic curvature. The self-accelerating branch of DGP model can explain the late time cosmic 
speed-up without recourse to DE or other components of energy mm- However, the self-accelerating DGP branch 
has ghost instabilities and it cannot realize phantom divide crossing by itself. To realize phantom divide crossing it 
is necessary to add at least a component of energy on the brane. On the other hand, the normal DGP branch cannot 
explain acceleration but it has the potential to realize a phantom-like phase by dynamical screening on the brane. 

In the present work we consider the HDE model in the framework of DGP braneword with GO cutoff, L = 
(aH + /3Ff 2 ) -1 / 2 , proposed in 0. Our work differs from 0 in that, they studied HDE model with GO cutoff in 
standard cosmology, while we investigate this model in the framework of DGP braneword and incorporate the effect 
of the extra dimension on the evolution of the cosmological parameters on the brane. The main difference between 
the HDE with GO cutoff in the framework of DGP braneword, with the one considered in standard cosmology 0 ], 
is that the equation of state parameter of the HDE with GO cutoff in standard cosmology is a constant 0 , namely 


w .d = -1 


2(a~l) 
3/3 ’ 


(3) 


however, as we shall see in DGP braneworld, due to the bulk effects, wd becomes a time variable parameter. Clearly, 
a time variable DE is more compatible with observations. In particular, the analysis of data from WMAP9 or Planck- 
2013 results on CMB anisotropy, BAO distance ratios from recent galaxy surveys, magnitude-redshift relations for 
distant SNe la from SNLS3 and Union2.1 samples indicate that the time varying DE gives a better fit than a 
cosmological constant [22], [0 . In addition, current data still slightly favor the quintom DE scenario with EoS across 
the cosmological constant boundary % = — 1 j22j. Furthermore, we are able to derive explicitly, the cosmological 
parameters on the brane as functions of redshift parameter and provide a profile of the cosmic evolution on the 
brane. Indeed, this is the advantages of the present model in compared to all other HDE models. As far as we know, 
this is the first model of HDE which leads to analytical solution, w = w(z) and q = q(z ). Then, we establish the 
correspondence between our model and scalar field models of DE and reconstruct the evolution of scalar fields and 
potentials. Finally, we investigate the stability of this model against perturbation in different cases. We find that for 
some range of the parameter spaces our model is stable which indicates the viability of this model for explanation of 
the late time acceleration. 


II. HDE IN DGP BRANEWORLD 

We consider a homogeneous and isotropic FRW universe on the brane which is described by the line element 


ds 2 = — df 2 + a 2 (t) 


dr 2 


1 — kr 2 


r 2 dH 2 


(4) 


where k = 0,1,—1 represent a flat, closed and open maximally symmetric space on the brane, respectively. The 
modified Friedmann equation in DGP braneworld model is given by [0 


H z 


k 


/ 9 

1 

3 M 2 

4r 2 


1 e \ 2 
r 2 2 r r ) ’ 


(5) 


where e = ±1 corresponds to the two branches of solutio ns [2(j| . and r c stands for crossover length scale between the 
small and large distances in DGP braneworld defined as (20j 


M pi _ g 5 


2 Mf 2G 4 


r, 


( 6 ) 
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We also assume that there is no energy exchange between the brane and the bulk and so the energy conservation 
equation holds on the brane, 


p + 3H(1 + w)p = 0. 


(7) 


For r c 1, the Friedmann equation in standard cosmology is recovered 


2 fc _ _p_ 
a 2 3M p Y 


( 8 ) 


Recent observations indicate that our universe is spatially flat. For a flat FRW universe on the brane, Eq. (0 reduces 
to 


H ~ r r H 3 Ml 


pi 


(9) 


Depending on the sign of e, there are two different branches for the DGP model. For e = +1 and in the absence of 
any kind of energy or matter field on the brane (p = 0), there is a de-Sitter solution for Eq. (0 with constant Hubble 
parameter 


H = 


a(t ) = aoe r c. 


( 10 ) 


Clearly, Eq. (1101) leads to an accelerating universe with constant equation of state parameter ujjj = — 1, exactly like the 
cosmological constant. However, there are some unsatisfactory problems with this solution. First of all, it suffers the 
well-known cosmological constant problems namely, the fine-tuning and the coincidence problems. Besides, it leads to 
a constant cod, while many cosmological evidences, especially the anal ysis of the type la supernova data indicates that 
the time varying DE gives a better fit than a cosmological constant [22h24| . Most of these data favor the evolution 
of the equation of state parameter and in particular it can have a transition from ujd > — 1 to ojd < — 1 at recent 
stag e. Although some evidence such as the galaxy cluster gas mass fraction data do not support the time-var ying wd 
[ 251] . an overwhelming flood of papers has appeared to understand the uijj = —1 crossing in the past decade [261 ]. In 
addition, to arrive at Eq. (HI one ignores all parts of energy on the brane including DE, dark matter and byronic 
matter, which is not a reasonable assumption. 

On the other hand, for e = — 1 and r c <C one can neglect the term H 2 in Eq. m to arrive at 


H 2 


36Mf' 


( 11 ) 


This is the Friedmann equation in spatially flat RS II braneworld [27|. Clearly, Eq. m does not have a self 
accelerating solution, so it implies the requirement of some kind of DE on the brane. 

In the present paper, we consider the HDE model in flat DGP braneworld with GO cutoff, which is defined as fl3| 


L={aH 2 + /3H)~ 1/2 . 


( 12 ) 


With this IR cutoff, the energy density © can be written 


Pd = 3M 2 (aff 2 + /3H), 


(13) 


where a and /3 are constants which should be constrained by the recent observational data and we have absorbed 
the constant c 2 in a and /?. Hereafter, we work in a unit in which M pl = 1. We also restrict our study to the current 
cosmological epoch, and hence we are not considering the contributions from matter and radiation by assuming that 
the dark energy pjj dominates, thus the Friedman equation becomes simpler. Substituting Eq. m into Friedmann 
equation 0 we can obtain the differential equation for the Hubble parameter as 


H 2 (l~a)- -H-/3H = 0. 
r c 

Solving this equation, the Hubble parameter is obtained as 


(14) 


r c (l — a) + ciee^= 


H(t) = 


(15) 
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where c\ is constant of integration. Since for r c 1, the effects of the extra dimension should be disappeared and 
the result of [Tdj . namely 


m = 


p i 


a — It 

must be restored, thus the constant c\ should be chosen as 

r c {a - 1) 


Cl = 


Substituting ci in Eq. m, we obtain 


H(t) = 


Taking the time derivative of Eq. m yields 


m = 


r c ( 1 — a)(l — ) 


eP r c 


Prl( 1 - a) (i _ e /fc)2' 


We can also solve Eq. (1181) to obtain the scale factor. We find 

a(t) = ao — 1 


ot-1 


(16) 


(17) 


(18) 


(19) 


( 20 ) 


Using the fact that 1 + z = ao/a, where 2 is the redsliift parameter, and combining Eqs. (1181) and (EU1) we find 
explicitly the Hubble parameter as a function of z, 


H{z) = 


l + (l + 2)“ 


rc(l-a) (i + z)^ 

£ r i + (l + 2)‘ 


r c (l - a) I 


( 21 ) 


This equation is valid in two cases. First, for e = +1 and a < 1, and second for e = — 1 and a > 1, and has no solution 
for a = 1. Inserting Eqs. m and m in m we get 


3 a + (1 — o)e' 3r <= 
PD r2(l-a) 2 (l_e^) 2 


Substituting Eqs. m and 03 in the time derivative of Eq. o. 

p D = 3(2aHH + j3H), 


( 22 ) 


(23) 


we arrive at 


Combining Eqs. 
function of time, 


Pd = 


3eef> r 


2(a + (1 — a)eP r c ) 


+ 1 — Q! 


(24) 


/3rj)(l — o) 2 (1 — ) 2 L (1 — ) 

03, 03 and (l24l) with 0, we find the equation of state parameter of HDE on the brane as a 


Wfl(f) = -1 - 


(l-o)c 


Pr c 


3/3 


: + (1 — o)e 


Pr c - 


(25) 


In order to see the evolution of u>d during the history of the universe, it is better to express it as a function of the 
redshift parameter 2 , 


u D (z) = -1 - 


1 —or (1 + a)(l + 2 )“^ +2 

3ft l-(i + (i + 2)~p~)(i + o(i + z)-^~y 


(26) 

























5 




FIG. 1: The evolution of u>d versus redshift parameter z for different model parameters. Left panel corresponds to a < 1 and 
(3 < 0, while the right panel shows the case a > 1 and /3 > 0. 


To check the limit of E7>1) in standard cosmology where r c 1, we note that, from Eq. COD and relation a/ao = 
(1 + z) _1 we have (a/ao) 2 ^ - = (1 + —» 0, asr c > 1. Therefore, ECU) reduces to 

2a- 1 


wd = — 1 + 


3 P 


(27) 


which is exactly the result obtained in jl3i |. It is worth mentioning that unlike in standard cosmology, where the 
equation of state parameter of HDE with GO cutoff becomes a constant 0 , in DGP braneworld scenario uid varies 
with time. Thus, it seems that the presence of the extra dimension brings rich physics. For a = 1 we have u>d = —1, 
similar to the cosmological constant. Let us consider the case with a < 1 and a > 1 separately. In the first case where 
a < 1, we find that the equation of state parameter can explain the acceleration of the universe provided p < 0. In 
the second case where a > 1, the accelerated expansion can be achieved for /? > 0. In both cases, our universe has 
a transition from deceleration to the acceleration phase around 0.4 < z < 1, compatible with the recent observations 
[28l - l30| . and mimics the cosmological constant at the late time. These results can be easily seen in Fig. (T] which shows 
the behaviour of ojd versus z for both cases. 

The first and second derivatives of the distance can be combined to obtain the deceleration parameter q. It was 
shown that the zero redshift value of go, is independent of space curvature, and can be obtained from the first and 
second derivatives of the coordinate distance [28| . It was argued that qo, which indicates whether the universe is 
accelerating at the current epoch, can be obtained directly from the supernova and radio galaxy data 28]. The 
deceleration parameter is given by 


Q = 



Using Eqs. m and m, the deceleration parameter is obtained as 


q = ~ 1 - 


1 — a 


l + (l + z) 



(28) 


(29) 


For a = 1 we have q = — 1. Although, the equation of state and the deceleration parameters do not depend 
explicitly on the crossover length scale r c which is the characterization of the DGP branworld, they depend on r c 
via the relation between the scale factor and the redshift parameter as in (EU1) . We have plotted the behavior of the 
deceleration parameter versus z in figure [2] Form this figures, we see that the transition from deceleration phase to 
the acceleration phase can be occurred around 0.4 < z < 1, which is consistent with the recent observations [281 ]. The 
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zero-redsliift value of the deceleration parameter is obtained as 

«=~ 1 + ^W- ( 3 °) 

Recent cosmological data from the combined sample of 192 supernovae and 30 radio galaxies, show that q 0 = —0.30 ± 
0.18, for a window function of width 0.4 in redshift and a transition from deceleration to the acceleration phase at 
redshift 2 = 0.78 ± 0.37 J28|. The zero-redshift value (1501) . for a = 0.7 and /3 = —0.22 leads to qo = —0.32, while for 
a = 1.2 and /3 = 0.15, we get go = —0.33, which is consistent with the observations [281 1 - 




FIG. 2: The evolution of the deceleration parameter q versus redshift parameter 2 for different model parameters. Left panel 
corresponds to a < 1 and /? < 0, while the right panel shows the case a > 1 and /3 > 0. 


III. CORRESPONDENCE WITH SCALAR FIELD MODELS 


The dynamical DE proposal is often realized by some scalar field mechanism which suggests that the energy form 
with negative pressure is provided by a scalar field evolving down a proper potential. Scalar fields naturally arise in 
particle physics including supersymmetric field theories and string/M theory. Therefore, scalar field is expected to 
reveal the dynamical mechanism and the nature of DE. However, although fundamental theories such as string/M 
theory do provide a number of possible candidates for scalar fields, they do not predict its potential, V((f>), uniquely. 
Consequently, it is meaningful to reconstruct the potential V (</>) from some DE models possessing some significant 
features of the quantum gravity theory, such as HDE models. Famous examples of scalar field DE_models include 


quintessence [3l|, K-essence [32j, tacliyon |33(, phantom [34J, ghost condensate 


quintom 37], and so forth. 


For a comprehensive review on scalar filed models of dark energy, see [38l |39{ . Generically, there are two points of 
view on the scalar field models of dynamical DE. One viewpoint regards the scalar field as a fundamental field of 
the nature. The nature of DE is, according to this viewpoint, completely attributed to some fundamental scalar field 
which is omnipresent in supersymmetric field theories and in string/M theory. The other viewpoint supports that 
the scalar field model is an effective description of an underlying theory of DE. If we regard the scalar field model as 
an effective description of such a theory, we should be capable of using the scalar field model to mimic the evolving 
behavior of the HDE and reconstructing the scalar field model according to the evolutionary behavior of HDE. 

In this section we implement a correspondence between HDE in DGP braneworld with GO cutoff and various 
scalar field models, by equating the equation of state parameters for these models with the obtained equation of state 
parameter of Eq. (EH). 
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A. Reconstructing holographic quintessence model 


Let us start with reconstructing the potential and dynamics of quintessence scaler field. The energy density and 
pressure of the quintessence scalar field are given by 


P4> = \<t> 2 + v {4>), (31) 

P4> = \^ % -V(<!>)■ (32) 

Thus the potential and the kinetic energy term can be written as 

Vffl = (33) 

4> 2 = (1 + (34) 


where ui$ = p^/' p$. In order to implement the correspondence between HDE and quintessence scaler field, we identify 
P<t> = Pd and = u>d- Using Eqs. (1251) . (1351) and (1551) . we find 


(1 + a)(l + z)—p~ + 2 


/3r 2 c {a-l) (1 + *) 


2(l-a ) 


(35) 


V(z) = 


2(a — l) 

3 + 


2 r 2 (1 - a) 2 

Eq. (1551) can also be rewritten 


(2a(l + z )'V + 2) [(1 + z) V + l] + I- ((1 + a)(l + z )V + 2 ) 


3/3 


(36) 


d<j> 

din a 


= = 


1 — a 


(1 + a)(l + z)—p~ + 2 


z^//3(a-l) 1 + (1 + 2;)" 


Using relation a = ao(l + z) 1 one can also obtain 


Integrating, yields 


d(f> 1 la— 1 ^/(l + <a)(l + z) p +2 

d2 e(l + z)V ^ l + (l + z) i ? i 


(37) 


(38) 


(j){z) 


dz f^T \!(1 + q)(l + z)^ +2 

e(l + 2)V /3 i + 


(39) 


Therefore, we have established holographic quintessence DE model and reconstructed the potential and the dynamics 
of scalar field according to evolutionary behaviour of HDE on the brane. Theoretically, one may expect to omit the 
redshift z from (156)) and (1391) to derive the potential as a function of the scalar field, namely, V(<j>). However, due to 
the complexity of the equations, the analytical form of the potential V = is hard to be derived. Nevertheless, 

we can plot the reconstructed potential as a function of (j> numerically. The reconstructed quintessence potential V ((/>) 
and the evolutionary form of the field are plotted in figure [5) where we have taken the zero redshift value of the scalar 
field equal to zero, namely, <j>(z = 0) = 0. Selected curves are plotted for a < 1 and the different model parameter 
/3 < 0. We can also plot the figures for the case a > 1 and /3 > 0, however, the behaviour is the same as the former 
case and for the economic reason we have not plotted here the latter case. From these figures we can see the dynamics 
of the potential as well as the scalar field explicitly. Figure [5] shows that the reconstructed quintessence potential is 
steeper in the early epoch and becomes flat near the present time. In other words, it mimics a cosmological constant 
at the the present time. 
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FIG. 3: Left panel shows the reconstruction of the quintessence potential V{<j>) for different parameters. In the right panel the 
behavior of the quintessence scalar field (p(z) for HDE model in DGP braneworld is illustrated for different parameters. 


B. Reconstructing holographic tachyon model 

The tachyon field is another candidate for DE. The equation of state parameter of the rolling tachyon smoothly 
interpolates between —1 and 0 (4flj |. Thus, tachyon can be realized as a suitable candidate for the inflation at high 
energy [idj as well as a source of DE depending on the form of the tachyon potential (42j. Choosing different self¬ 
interacting potentials in the tachyon field model lead to different consequences for the resulting DE model. Due to all 
the above reasons, the reconstruction of tachyon potential V(<j>) is of great importance. The correspondence between 
tachyon field and various DE scenarios such as HDE [43| and agegraphic DE [44] has been already established. The 
study has also been generalized to the entropy corrected holographic and agegraphic DE models j45j . 

The tachyon condensates in a class of string theories and can be described by an effective scalar field with a 
Lagrangian of the form by [ 1 (> 


L = - g^d^dv^), (40) 

where V (</>) is the tachyon potential. The corresponding energy momentum tensor for the tachyon field can be written 
in a perfect fluid form 


= {p + p^pU,, - pgfj, v 

where p and p are, respectively, the energy density and pressure of the tachyon, and the velocity u^ is 

The energy density and pressure of tachyon field are given by 


l-i 


p = -t 0 ° = 

Thus the equation of state parameter of tachyon field is given by 

P 12 , 

U>T = — =</> — 1 . 

P 


(41) 


(42) 


(43) 

(44) 


( 45 ) 
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To establish the correspondence between HDE and tachyon field, we equate ojd with u>t- Combining Eqs. m, © 
and (TH1) . we find 


V(z) 

4> 2 {z) 


3 

1 

(l + a)(l + js)V +2 

(1 + a(l + z)~r L )(1 + (1 + z)~^) 

r c(l - a) 2 \ 

3/3 

_(! + (! + z)^~)( 1 + a(l + z)^~) 

. 2(l-a) 

(1 +Z)—P— 


, (46) 


a — 1 

~W 


(1 + a)(l + z)~^~ + 2 
_(! + (! + 0) _ ^ _ )(1 + a{l + z) 



(47) 


We can further rewrite {HD as 


d< t> H - x 6- - r °( 1 + z ^ .) 

/(a-l) 3 

(l + a)(l + z)V +2 

d(lno) * e Y 

3/3 

_(! + (! + z) _ ?') 3 (l + a(l + z)~^)_ 


Using the fact that 


d(j) d(f> da 
dz da dz 


-(1 + *) 


-i d(j) 


d(lna) ’ 


(48) 


(49) 


we have 


d(f> 

r c (l + z)~^ 

/(a - l) 3 

(l + a)(l + 2)V +2 

dz 

e(l + z) 

V 3/3 

_(! + (! + z)~^) 3(1 + a(l + z)~ir L )_ 


(50) 


Integrating and setting the constant of integration equal to zero by assuming <j>{z = 0) = 0, we obtain 
<t>{z) = I ' — - —\! _ y - dz. 


f r c ( 1 + z)~^ 

/(a-l) 3 

(1 + 

a)(l + z)~^ + 2 

/ e(l + z) 

V 3/9 

.(! + (!+*) 

V)3(l + a (l + .)V)_ 


(51) 


In this way we connect the HDE in DGP braneworld with a tachyon field, and reconstruct the potential and the 




FIG. 4: Left panel shows the reconstruction of the tachyon potential V(tj>) for different parameters. In the right panel the 
behavior of the tachyon field <f>(z) for HDE model in DGP braneworld is illustrated for different parameters. 


dynamics of the tachyon field which describe tachyon cosmology. The reconstructed tachyon field and the evolution 
of the tachyon potential are plotted in figure [I] Again, we see that tachyon potential is steeper in the early epoch 
and becomes flat near today. Thus, the universe enters a de-Sitter phase at the late time and the constant potential 
plays the role of cosmological constant. 
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IV. STABILITY OF THE MODEL 

In this section we would like to investigate the stability of the HDE in DGP braneworld against perturbation. It 
is expected that any viable DE model should result a stable DE dominated universe. A simple, but not complete, 
approach to check the stability of a proposed DE model is to study the behavior of the square sound speed (v 2 = dp/dp) 
[47]. It was argued that the sign of v 2 plays a crucial role in determining the stability of the background evolution. 
If v 2 < 0, it implies that the perturbation of the background energy density is not an oscillatory function and may 
grow or decay with time, and so we have the classical instability of a given perturbation. On the other hand, the 
positivity of v 2 > 0 indicates that the perturbation in the energy density, propagates in the environment and so we 
expect a stable universe against perturbations. It is important to note that the positivity of v 2 is necessary but is 
not enough to conclude that the model is stable. Indeed, the negativity of it shows a sign of instability in the model. 
The behavior of the squared sound speed for HDE Idsl , agegraphic DE |ji) and the ghost DE model [13, HU were 
investigated. It was found that all these models [48j-l5l| are instable against background perturbations and so cannot 
lead to a stable DE dominated universe. 

In the linear perturbation regime, the perturbed energy density of the background can be written as 

p{t,x) = p{t) + Sp(t,x), (52) 

where p(t) is unperturbed background energy density. The energy conservation equation (V ll T IJ,v = 0) yields [47| 

Sp = v 2 S7 2 Sp{t,x). (53) 

For v 2 > 0, Eq. (l53l) becomes a regular wave equation whose solution is given by Sp = Spoe ~ lulot+lk ' x , which indicates a 
propagation mode for the density perturbations. For v 2 < 0, the perturbation becomes an irregular wave equation and 
in this case the frequency of the oscillations are pure imaginary and the density perturbation will grow with time as 
Sp = Spoe ut+lk ' x . Hence the negative squared speed shows an exponentially growing mode for a density perturbation. 
Thus the growing perturbation with time indicates a possible emergency of instabilities in the background. The 




FIG. 5: The behaviour of v 2 for HDE model in DGP braneworld for different model parameters. Left panel corresponds to 
a < 1 and /3 < 0, while the right panel shows the case a > 1 and /3 > 0. 


squared speed of sound can be written as 


v 


2 

S 


dp p 
dp P 


(54) 


where 


p = u: D p + u D p. 


( 55 ) 
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Taking the time derivative of (l26l) and substituting the result in (l54|) . after some calculation, we arrive at 


v 2 = 


a — 1 

~w~ 


4 + (1 + a)(l + z)~ 


2 +(1 + a)(l + zY 


- 1 . 


(56) 


The evolution of v 2 versus redshift parameter are plotted in figure [5] From these figures we see that v 2 can be positive 
provided the parameters of the model are chosen suitably. For example, for a < 1, /3 < 0, and also for a > 1, /3 > 0 
the square of sound speed is always positive during the history of the universe, and so in these cases the stable DE 
dominated universe can be achieved. 


V. CONCLUDING REMARKS 

In this paper, we have investigated the holographic model of DE in the framework of DGP braneworld. We have 
chosen the GO IR cutoff in the form, L = (aH + /37? 2 )” 1 / 2 , with two free parameters a and /?, which should be 
constrained by comparing with observations. This cutoff is the generalization of the well-known Ricci scalar cutoff 
in flat universe. We have restricted our study to the current cosmological epoch, and so we have not considered the 
contributions from matter and radiation by assuming that the dark energy pjj dominates, thus the Friedman equation 
becomes simpler. The main difference between the HDE with GO cutoff studied in this work in the framework of 
DGP braneword, with the one considered in standard cosmology fl3| . is that the equation of state parameter of 
the HDE with GO cutoff in standard cosmology is a constant [13|, however, in DGP braneworld, due to the bulk 
effects, wd becomes a time variable parameter. The application of the HDE with GO cutoff, in DGP braneworld 
allows us to solve the Friedmann equation and derive the Hubble parameter, H(z), as well as the scale factor, a(t ), 
analytically. We also obtained the equation of state and the deceleration parameters of HDE as the functions of 
the redshift parameter z and plotted the evolutionary behaviour of these parameters against 0 . The cosmological 
quantities depend on the two model parameter a and /?. To show the viability of the model, we studied the zero- 
redshift values of the cosmological parameters by chosen the suitable values for a and /3. For example, we found that 
for a = 0.7 and /? = —0.21, —0.22, —0.23, the transition from deceleration phase to the acceleration phase occurred 
around 0.4 < z < 1, which is consistent with recent observations (28|. The zero-redshift value of the deceleration 
parameter was obtained q 0 = —0.32 for a = 0.7, /3 = —0.22, and qo = —0.33 for a = 1.2, /3 = 0.15 which is again 
compatible with the cosmological data [28|. Besides, for a = 1 this model mimics a cosmological constant with 
u>d = —1, and q = —1, independent of the redshift parameter z. For a > 1, (3 > 0 and a < 1, (3 < 0 the equation of 
state parameter, however, is always larger than —1, and the Universe enters a de-Sitter phase at the late time. 

We have also established a connection between the quintessence/tachyon scalar held and the HDE in DGP 
braneworld. As a result, we reconstructed the corresponding potentials of the scalar held, V((/)), and the dynam¬ 
ics of the scalar helds as a function of redshift parameter, (f> = <(>( 2 ), according to the evolutionary behavior of the 
HDE model. Finally, we studied the stability of the presented model by studying the evolution of the squared sound 
speed v 2 whose sign determines the sound stability of the model. Interestingly enough, we found that v 2 > 0 provided 
the parameters of the model are chosen suitably. For example, for a = 0.7, /3 = —0.1, and a = 1.2, (3 = 0.06 the 
squared sound speed is always positive during the history of the universe, and so in this case the stable DE dominated 
universe can be achieved. This is in contrast to HDE in standard cosmology, which is instable against background 
perturbations and so cannot lead to a stable DE dominated universe j48| . This implies that the presence of the extra 
dimension in HDE model, can bring rich physics and in particular it has an important effect on the stability of the 
HDE model. This issue deserves further investigations. 
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